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Abstract 

The techniques for polynomial interpolation and Gaussian quadrature are gen- 
eralized to matrix- valued functions. It is shown how the zeros and rootvectors of 
matrix orthonormal polynomials can be used to get a quadrature formula with the 
highest degree of precision. 

1 Introduction 

The aim of this paper is to construct quadrature formulas, using orthogonal matrix poly- 
nomials, to approximate matrix integrals. We will give an expression for the quadrature 
coefficients and show that the constructed formula has the highest possible degree of 
accuracy and converges to the exact value of the matrix integral. All these ideas are 
generalizations of the classical Gaussian quadrature rules for the scalar case. 

In Section 2 we will give a survey of definitions and properties of matrix polynomials. 
These can be found in the book on matrix polynomials by I. Gohberg, P. Lancaster and L. 



Rodman [|1(J and in the survey on orthogonal matrix polynomials by L. Rodman []12}| . In 
Section 3 we introduce orthogonal matrix polynomials on the real line and discuss some 
properties which we will need in the next sections. These orthogonal matrix polynomials 
have been considered earlier by Delsarte, Genin and Kamp [|3] and Geronimo ||, J/J. As 
in the scalar case, the theory of approximate integration uses results from the theory 
of interpolation. In Section 4 we will discuss polynomial interpolation and particularly 
the interpolation problem of Lagrange. We will give an expression for the Lagrange 
interpolation polynomial in the general case and then apply this result to the case of 
a Jordan pair (A, J) of the orthonormal matrix polynomial P n {x). These results are 
also generalizations of the known results for the scalar case. The interpolation problem 
has also been treated in ||, ||, |3J] and but their approach is an algebraic one. We 
have restricted our attention to the interpolation formula we needed for the construction 
of Gaussian quadrature rules. The Gaussian quadrature formula is then constructed in 
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Section 5, where we give a formula for the quadrature coefficients and show that the 
quadrature rule converges under appropriate conditions. 



2 Matrix polynomials 

If Aq, Ax, . . . , A n are elements of MP xp and A n 7^ 0, then we call 

P(x) = A n x n + A n _ x x n ~ l + . . . + A x x + A , 

a matrix polynomial of degree n. This matrix polynomial is monic when A n = I, the 
identity matrix. A point Xq is a zero of P{x), if det P(xq) = 0. Note that if the lead- 
ing coefficient of P(x) is non-singular, detP(x) is a polynomial of degree np. Another 
important notion associated with matrix polynomials are Jordan chains. A sequence of 
p-dimensional column vectors v ,vi, . . . ,Vk is called a right Jordan chain of length k + 1 
of a monic matrix polynomial P(x) corresponding to Xq, if v 7^ and 

1 1 

^-P w (x )^ = 0, Z = 0,l,...,fc. 

2 = l - 

The initial vector vq 7^ is called a rootvector of P(x) corresponding to xq. Note that in 



flTOjj and [|T^] the zeros are called eigenvalues and the rootvectors are called eigenvectors. 
p( l \x) is the ith. derivative of P(x) with respect to x and this means that we take the 
ith derivative of every element of P(x) with respect to x. 

In analogy with the definition of a right Jordan chain, we call a sequence of p- 
dimensional row vectors Wq, Wi, . . . , a left Jordan chain of length k + 1 of a monic 
matrix polynomial P(x) corresponding to Xq, if Vq 7^ and 

1 1 

^-^pW(x ) = 0, Z = 0,l,...,fc. 
i=o % - 

Jordan chains are not unique: a matrix polynomial can have different Jordan chains 
of various lengths and different rootvectors, corresponding to the same zero. In what 
follows we formulate the definitions and properties for right Jordan chains, unless explicitly 
mentioned. 

A convenient way of writing a Jordan chain is given by the following property. 



Proposition 2.1 ( |10[| , p. 27) The vectors Vq, v 1, . . . , Vk form a right Jordan chain of 



the monic matrix polynomial P(x) = Ix n + A n _\x n 1 + . . . + A\x + Aq corresponding to 
xo if and only if v 7^ and 

Ao Jq + A^iXqJq 1 + . . . + AiXqJo + AoXq = 0, 

where Xq = v ... Vk ) is a p x (k + 1) matrix and Jq is a Jordan block of size 
(k + 1) x (k + 1) with xq on the main diagonal. 

Observe that the equations 



y £~L®(x o )v l - i = i 



0,1,..., k, 
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where L(x) = L„x r ' 



L n -ix n 1 + . . . + L\x + L , with Li e 



ppx<j 



(i=l,2,. . . ,n), can always 



be written as 

LuXqJq + L n _iX Jq + . . . + LiXqJ + L X = 0, 

where X = ( vq ... % J is a p x (H 1) matrix and Jo is a Jordan block of size 
(k + 1) x (A; + 1) with xo on the main diagonal. 

In the following definitions and properties we restrict ourselves to monic p x p matrix 
polynomials, but most of the theory can also be given in the context of regular matrix 
polynomials. These are matrix polynomials which satisfy det P(x) ^ 0. 

Now we introduce the notion of a canonical set of Jordan chains. Let 



.(*) .,(0 



.(*) 



1,2, 



be a set of Jordan chains of a p x p monic matrix polynomial P(x) corresponding to the 



zero Xi. Then we call the set canonical if the rootvectors Uioi^oj 



, v® are linearly 

independent and J^jLi^f — m ii where is the multiplicity of X{ as zero of P(x). Such 
a canonical set of Jordan chains is not unique, but the number of chains and their length 
depend only upon P(x) and Xi and do not depend on the choice of canonical set. A 
canonical set of Jordan chains can be associated with a pair of matrices (Xi, Jj), which is 
called the Jordan pair of P(x) corresponding to Xj and defined as follows : 



Xi 



(v {i) 

\ v l,0 



(') 



(0 



(0 
•r ' 



(0 



a p x rrii dimensional matrix 



and 



Ji 



diag(J itll J ij2 , 

I Xi 1 



Ji 



l, Si 



a rrii x rrii dimensional matrix with 



V 



1 

Xi ) 



a /ij x dimensional matrix. 



A pair of matrices (X, J) where X is a p x np dimensional matrix and J a np x np 
dimensional Jordan matrix is called a Jordan pair for the monic matrix polynomial P(x) 
if 

X = (Xi X 2 . . . X k ) and J = diag(J\ J 2 ■ ■ ■ Jk), 

where (Xi, Ji) is a Jordan pair of P(x) corresponding to Xi and k is the number of different 
zeros of P(x). Jordan pairs have the following important property : 

Proposition 2.2 ( ||10|| , p. 45) Let (X, J) be a pair of matrices where X is of size px np 
and J is a Jordan matrix of size np x np. Then (X, J) is a Jordan pair of the monic 
matrix polynomial P(x) = Ix n + A n ^\x n ~ x + . . . + A\x + A if and only if 

(1) 



( 



ool(XJ l )Zl 



X 
XJ 
XJ 2 



\ xj n - L J 



is a non-singular np x np matrix, 
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(2) XJ n + A n _ x XJ n ~ l + ... + A 1 XJ + A X = 0. 



The requirement that J is Jordan is not essential. We call a pair of matrices (X, T) 
where X is of size p x np and T is a np x np dimensional matrix, a standard pair for the 
monic matrix polynomial P(x) = Ix n + A n _\x n ~ x + . . . + A\x + A if 



(1) col(XT l )™ =0 is a non-singular np x matrix, 

(2) XT n + A^iir- 1 + . . . + AiXT + A X = 0. 

This means that every Jordan pair is a standard pair and every standard pair (X, T) for 
which T is a Jordan matrix is a Jordan pair. 

With every standard pair (X, T) of a monic matrix polynomial P(x) we can associate 
a third matrix Y of size np x p, with 



y 



x 

XJ 
XJ 2 



\ _1 / o \ 





XJ 



n-l 





V 1 J 



The triple (X,T,Y) is called a standard triple for P(x) and if T = J, a Jordan matrix, 
then (X, J, y) is called a Jordan triple. One can proof that (X', Ci, Y'), where 



X' 



( / ... ) , d 



V 








-A 



I 








/ 



-A 2 










/ 

-A n -! J 



and Y' 



( \ 




V i ) 



is a standard triple of P(x) = Ix n + A„_ 1 x™ _1 + . . . + A x x + A . Two standard triples 
(X, T, y) and (X' ,T' ,Y') are similar if there exists an invertible np x np matrix S such 
that 

X' = XS ; T' = 5 _1 T5 and Y' = S^Y. 



This matrix S is uniquely defined by, 



S = (col{XT l )"J) 'col{X'T 



Consider a Jordan triple (X, J, y), then we already know that the columns of X, when 
decomposed into blocks consistently with the decomposition of J into blocks, form right 
Jordan chains for P(x). In analogy with this, we can give a similar meaning to the rows 
of y. Indeed, the rows of Y, partitioned into blocks consistently with the decomposition 
of J into blocks and taken in reverse order, form left Jordan chains for P(x). The notion 
of standard triples is important for the following representation theorem: 

Theorem 2.3 (IffOfl, p. 58) Let P(x) = Ix n + A n „ x x n ~ l + ■ ■ ■ + A x x + A be a monic 
matrix polynomial of degree n with standard triple (X,T,Y). Then P(x) admits the 
following representations : 
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(1) P{x) = x n I — XT n (V\ + V2X + . . . + V n x n 1 ), where Vi are npx p matrices such that 
(V, V 2 ... V n ) = {col{XT l )^f \ 

(2) P(x) = x n I — (W\ + W2X + . . . + W n x n ~ r )T n Y, where Wi are p x np matrices such 
that col(Wi)" =1 = (y TY T 2 Y ... T n ' x Y ) _i . 

Note that those forms are independent of the choice of the standard triple. 

Finally we have some properties about the divisibility of matrix polynomials. We say 
that the matrix polynomials Q(x) and R(x) are the right quotient and the right remainder, 
respectively, of P{x) on division by D(x) if 

P(x) = Q(x)D(x) +R(x) 

and if the degree of R(x) is less than that of D(x). The right division of matrix polynomials 
of the same order is always possible and unique, provided the divisor is a polynomial with 
non- singular leading coefficient (see |5|, p. 78]). In a similar manner we can define the left 
division : P(x) = D(x)Q(x) + R(x). 

Proposition 2.4 Let P(x) = A m x m + A m _\x m ~ x + . . . + Aix + A be a matrix polynomial 
of degree m and let D{x) = B n x n + B n _\x n ~ x + . . . + B\x + B be a matrix polynomial of 
degree n, with a non-singular leading coefficient B n and with Jordan pair (X, J). Then 
D(x) is a right divisor of P(x) if and only if 

A m XJ m + A^XJ" 1 - 1 + ... + A X XJ + A Q X = 0. 

Proof 

First of all we denote that every Jordan chain Vq, Vi, . . . , Vf. of D(x) and corresponding to 
Xq is also a Jordan chain of Q(x)D(x), corresponding to the same zero : 

EkQ(x)D(x)}<§ =x v^ = j:t l U i .)Q it) (x )D( i - t \x )v^ 
i=0 % - i=0 t=o 11 \ 1 ) 

= E^ (i) (^o)E 7 -^ T ^ ) (^o)^ 

t=0 s=0 S ' 

= for Z = 0,l,...,fc. 

Suppose D(x) is a right divisor of P{x), this means P(x) = Q{x)D(x). But then every 
Jordan chain of D(x) is also a Jordan chain of P(x). In particular we have 

A m XJ m + A m _iXJ m ~ l + ... + A X XJ + A X = 0. 

Suppose now we know that the above mentioned equation holds. This means that every 
Jordan chain of the Jordan pair of D(x) is also a Jordan chain of P(x). Moreover, since 
the leading coefficient of D(x) is non-singular, the right division is possible and unique, 
P(x) = Q(x)D(x) + R(x). So, every Jordan chain of the Jordan pair of D(x) is also a 
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Jordan chain of R(x). But this implies that the matrix polynomial R(x) of degree < n — 1 
satisfies 

( X ^ 

/ x XJ 

y R R\ ... R n -\ 



0. 



V xj n ~ x j 



Since (X, J) is a Jordan pair, the np x rap dimensional matrix col(X J 1 )™ = q is non-singular 
and thus R(x) = 0. ■ 

Note that this proposition was also proved in [^, Thm. 2.1 on p. 333]. 



3 Orthogonal matrix polynomials on the real line 

A symmetric p x p matrix- valued function W(x), integrable over [a, b] is called a weight 
matrix function if W(x) > and det W(x) ^ holds almost everywhere (see |§). The 
notation X < Y for symmetric matrices means that Y — X is positive semidefinite. The 
matrix integral 

/ F(x)W(x)G(x) T dx 

J a 

where F(x) and G(x) are continuous matrix-valued functions, is defined in a natural way. 
The (i,j)th element is given by a sum of integrals : 

EE [ b F(x). s W(x) Sjt G(x) T ttj dx. 
s= it=i Ja 

Let R pxp [x] be the set of polynomials in a real variable x and whose coefficients are p x p 
matrices with real entries. If P(x) and Q(x) are elements of this set, then we define a 
matrical inner product on R pxp [x] as follows: 

(P(x),Q{x)) L = [ P{x)W{x)Q{x) T dx. 

J a 

This matricial inner product has some properties which we will recognize as generalizations 
of the properties of the scalar inner product. 

Proposition 3.1 1. (P,Q) L = (Q,P) L T where P,Q eW xp [x}. 

2. (CxPx + C 2 P 2 , Q) L = d(Pi, Q) L + C 2 (P 2 , Q) L where C u C 2 G R pxp and P u P 2 , Q G 
« pxp [x}. 

3. (xP,Q) L = (P,xQ) L where P,Q G R pxp [x}. 

4- Let P(x) G M pxp [x] ; then (P,P) L is positive semidefinite and even positive definite 
if det P(x) ^ 0. 

5. Let P G W xp [x], then (P, P) L = if and only ifP = 0. 

These properties are easily proved by means of straightforward computation. 
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A generalization of the Gram-Schmidt orthonormalisation procedure for the set 
{I,xI,x 2 I, . . .} with respect to the matricial inner product (., .) L will give a set of or- 
thonormal matrix polynomials {Pn{%)}™=o which satisfy 

{ P n (x) W(x) P m {x) T dx = S n;m I. 

J a 

Moreover, P n (x) is a matrix polynomial of degree n, with a non-singular leading coefficient 
and is defined upon a multiplication on the left by an orthogonal matrix. 

As in the scalar case, these orthonormal matrix polynomials are orthogonal to every 
matrix polynomial of lower degree and they satisfy a three-term recurrence relation. 

xP n (x) = D n+1 P n+1 (x) + E n P n (x) + D 7 P n ^(x) } n > 0, 

P-i{x) = and P (x) = I, 

where D n is a positive definite matrix and E n is a symmetric matrix. The orthonormal 
polynomials are defined only up to a left orthogonal factor and it is convenient to choose 
this factor in such a way that the recurrence coefficients D n are symmetric. We assumed, 
without loss of generality that f% W(x) dx = I. Furthermore we have the Christoffel- 
Darboux formula : 

n 

P n {y) T D n+1 P n+1 {x) - P n+1 {y) T D n+1 P n {x) = {x-y) £ P^y) 7 ' P^x) . 

i=0 

(see 0). If we take x = y , we get 

P n (x) T D n+1 P n+1 (x) = P n+1 (x) T D n+1 P n (x) 

so that P n {x) T D n+ iP n+ i(x) is a symmetric matrix. By means of straightforward compu- 
tation we get the following equation : 

n 

J2Pi(x) T Pi(x) = PZ +1 (x)'D n+1 P n (x) - pT(x)'D n+1 P n+1 (x). 

i=0 

The matrix 

K n {x 1 y) = Y J P l {y) T P i {x) 

i=0 

is a positive definite matrix and we call it the reproducing kernel because of the following 
property. 

Proposition 3.2 Let II m (x) be a matrix polynomial of degree m <n, then 

(U m (x),K n (x,y)) L = U m (y) 

Proof 

If we write U m (x) in terms of the orthonormal matrix polynomials Pq{x), . . . , P m (x) : 

m 

n m (x) = Y, A iPi(x), 

i=0 



7 



we have for m < n 



(Il m (x),K n (x,y)) 1 



i=0 j=0 
m 

T, A i p i(y) 

i=0 

n m (y). ■ 



In the scalar case all the zeros of an orthonormal polynomial are simple. This is not 
the case for orthonormal matrix polynomials, but nevertheless we can proof a similar 
property. 

Proposition 3.3 The zeros of the orthonormal matrix polynomial P n (x) have a multi- 
plicity < p, where p is the size of the matrices. 

Proof 

Let Xq be a zero of P n (x) with multiplicity m > p. Consider a canonical set of right 
Jordan chains corresponding to x 

Vi,o,Vi,i,---,Vi,ixi-i, i = l,2,...,s. 

This means that the p-dimensional column vectors i^o, t>2,o> ■ ■ ■ ■> v s,o are linearly indepen- 
dent and J2t=i Hi — m (Section 1). Since m > p, there has to be a Jordan chain of length 
> 1. Suppose Vq and v\ are the two leading vectors of this chain, then they satisfy 

^(^o^o = and v ^ 0, 
Pn(x )v + P n (x )vi = 0. 

Using these equations, we get 

Vq K^^Xq, X )V = Vq P„ (X )' DnP^^Vo - Vq P^ixo)' D n P n (x )v 

= -v[P^(xo)D n P n - 1 (xo)v 
= -vJP^ 1 (x )D n P n (x )v 
= 0. 

But K n ^i(x ,x ) is a symmetric and positive definite matrix and vq 7^ 0. So the multi- 
plicity of xq as zero of P n (x) has to be < p. ■ 

Corollary. 

In this proof we showed that the length of a Jordan chain of P n (x) cannot be greater 
than 1. Thus a canonical set of right Jordan chains of P n (x) corresponding with a zero 
Xq consist of m linearly independent, non-zero p-dimensional column vectors, where m is 
the multiplicity of xq as zero of P n (x). 
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4 Polynomial interpolation 

4.1 Polynomial interpolation in general 

Consider apxp matrix-valued function F(x) and k different points x±, £ 2 , ■ ■ ■ , x^ with mul- 
tiplicity resp. mi, m 2 , . . . , m& where X^Li m i = np. With every point Xi, i — 1, 2, . . . , k, 
we associate a set of p-dimensional column vectors 

(i) (i) (i) (i) (i) (i) (i) (i) (i) 

^1,0)^1,1) • • • > U , „(0 ,^2,0)^2,1) •• (0 1' ••• ) U Si,0) U Si,l) ••• „(0 1 

where SjLi A*^ = and v^q, v^q, ■ ■ ■ , v^ Q are non-zero, linearly independent vectors. 
When we put these vectors in a p x np dimensional matrix, we get 

X = (Xi X 2 . . . X k ) where X t = („£> . . . v^ . . . „» ). 

The square x np matrix J is given by 

J = diag(J u J 2 , . . . , J fc ) where J< = diag(J i:1 , J ij2 , . . . , J i>ai ) 

/ ^ 1 \ 



and Jjj 



f 



Here Jj is a mj x m,j matrix and Jj j a ^ x /i^ matrix. We require that the matrix 

coI(XJ 1 )™ = q is non-singular. Note that this is a necessary and sufficient condition in 
order that (X, J) is a Jordan pair of a monic matrix polynomial of degree n (Section 1). 
We will construct a matrix polynomial P(x) = A n _ix n ~ 1 + . . . + A x x + A , which satisfies 



1=0 



1=0 



and this for 



q = 0, 1, . . . ,/4 l) - 1, £ = 1,2,...,^ and i = 1, 2, . . . , 



The coefficients of P(x) can be determined by solving a linear system of np 2 unknows and 
np 2 equations, because Y$=i Z)t=i = n P- 

Theorem 4.1 The general polynomial interpolation problem has a unique solution. 
Proof 

Let P\{x) and P2(x) be two solutions of the polynomial interpolation problem. Then 
Pi(x) — P 2 (:r) is a matrix polynomial of degree < n — 1 which satisfies 



E*Fi (O (*0-^ ,; (*i)] 



(0, 



(») 



0, 
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where 



q — 0, 1, ... , — 1, t = 1, 2, . . . , Si and i = 1,2, . . . , k. 



From this we get that the vectors 

(i) (i) ({) 



t — 1, 2, . . . , Si, i — 1, 2, . . . , k 



form a Jordan chain for the matrix polynomial P± (x) — Pi (x) = B n _ix n 1 + . . . + E>ix + B , 
corresponding to Xj. But this implies that 



B n -iX it J™ t 1 + B n _iX it J™ t z + . . . + BiX i:t J it + B X i:t — 0, 



jn-2 



where t = 1,2, ... ,Si and % = 1,2, ... ,k. The p x //^ dimensional matrix X iyt and the 

{ Xi 1 \ 



fj^ x /if' 1 dimensional matrix J it are given by 



X 



■i.t 



(0 (») 



and J, 



Xi 1 

Xi J 



Hence we get 
( B B 1 



B n -i 



( X ht \ 

XitJi,t 



0, t — 1,2, Si, % — 1,2, . . . ,k, 



or 



( B B\ ... -B n _i ) 



/ X \ 
XJ 



= 0. 



Since colfXJ')^ 1 is a non-singular matrix, we have 

Bi = for i = 1,2, 1. 
And this implies that P\(x) = Pi{x). ■ 
Corollary 

If F(x) is a matrix polynomial of degree < n — 1, then = F(x). 

The interpolation problem can be formulated in terms of vector polynomials instead 
of the sets of p-dimensional vectors. Define the p-dimensional vectors 



1=0 



and the p-dimensional vector polynomials 



Vt\x) = U«*o + U S ( X - X i) + V tM X ~ X if + ■■■+ , ( X ~ X i) 



,(*) ..,(*)/ 
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(i) / \ (i) | (i) / \ , (i) / 

Z t l X J — ^,0 ~+~ Z t,l\ X ~ X i) i Z t,2\ X 



Xi) 2 + . . . + 2 (<) w (X - Xi)^" 1 
i,/i t 1 



for g = 0, 1, . . . , /4 — 1, t — 1, 2, . . . , and z = 1, 2, . . . , k. In this case we are looking for 
a matrix polynomial P(x) which satisfies 

/ Qa-1 



— P(x)vt(x 



dx c 



a-l 



r 4 i} (x) 



dx c 



for a = 1, 2, ... , /4 , t = 1, 2, . . . , s, and i = 1, 2, . . . , k. In addition we know v[ l \xi) are 
non-zero and linearly independent vectors. 

If we set Si = 1, i = 1, 2, . . . , fc, we have /c different points xi, x 2 , . . . , x*,, vector poly- 
nomials «i (x), (x), . . . , Vi(x) for which Vi\xi) ^ 0, z = 1,2, ...,k and vector 
polynomials z[ (x),zj (x),...,z[ (x) and we want to find a matrix polynomial P(x) 
which satisfies 



& 



a-l 



— P( x )vf{x) 



dx' 







ia— i 



<9x 



a-l 



zP(x) 



, a = 1,2, . . . ,rrii, 



1,2, ... , k. 



This is the interpolation problem treated in |TJ. 

We will put /4 = 1 instead of Si — 1 in the interpolation problem that we will discuss 
and which we need for the construction of the Gaussian quadrature rule. 



4.2 The interpolation problem of Lagrange 

If all [ij = 1, we get the interpolation problem of Lagrange. This means that we have 
k different points x%, x 2 , . . . ,Xk with multiplicity resp. mi, m 2 , ■ ■ ■ , and Ya=i m i — n P- 
Every point Xj, i — 1, 2, . . . , k, is associated with a set of non-zero, linearly independent 
p-dimensional column vectors v it i,Vi >2 , 

X = (V, i ... vi mi ... v k ,i ■ ■ • Vk, mt .) and J = diag(xi, 



, Vi iTrH . Let X and J be defined as follows 



Xi, X2, 



,x 2 , 



Xfc , 



> XfcJ , 



mi 



Hi 2 



Vk,l ■ ■ ■ v k,m k 

then we require that the matrix col(XJ l )? Q is non-singular. 

We will construct the interpolation matrix polynomial P(x) of degree n — 1 such that 

P{xi)Vij = F(x i )v i j, j = l,...,m u i = l,...,k. 

Special case. 

If p = 1, then all the multiplicities have to be equal to 1. This means we have n different 
points xi, x 2 , . . . , x n and every point Xj is associated with a non-zero number v i: i = 
1,2, ... ,k. The matrix 



X 
XJ 
XJ 2 



\ XJ 



n-1 



ViXi 
V\x\ 



V-2 
V2X2 

v 2 x\ 



V n X T[ 
2 

Vn.Xr 



J 


\ v x x\ 


1 „, ~n-l 
V2X2 


• • • v n x™ 














( 1 


1 


1 \ 









. 


. 


\ 




Xi 


X 2 


Xn 









. 


. 






,7-* 2 


™2 

x 2 • • • 












^3 • 


. 






K xT l 


X 2 ••• 


x™- 1 y 




V 





. 


• v n 


) 
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is non-singular because of our choice of the points xi, x 2 , ■ ■ ■ , x n and the non zero-numbers 



Vi,V 2 



1 *->- A , . . . , u n . 



The interpolation polynomial satisfies 



p{xi)vi = f{xi)vi, i = l,2,...,n. 



But this is equivalent with 

P(xi) = f(xi), i = l,2,...,n, 

because i>j 7^ 0, for i — 1,2, ...,n. So we find the polynomial interpolation problem of 
Lagrange, for real functions / : M — » R. 

If rrii = p ■ di, i = 1, 2, . . . , k and col(Xi J\)i 1 = q are non-singular matrices, then (Xj, Jj) 
is the Jordan pair of a monic matrix polynomial Li(x) of degree d{. Because di > 1, 
Ri(x) = F(xi) will be a matrix polynomial of degree < ci, and this for i = 1,2, ... ,k. 
Now we get the following interpolation problem: given Li(x), L 2 (x), . . . , L^x), monic ma- 
trix polynomials of degree d\, d 2 , . . . , d k respectively and Ri(x), R 2 (x), . . . , Rk(x) matrix 
polynomials of degree < di respectively. Find a matrix polynomial P(x) such that 



P(x) = Si(x)LAx) + RAx) 



1,2, 



k 



for some matrix polynomials S±(x), S 2 (x), ... , Sk(x). This interpolation problem is treated 
in 1, i and 

Theorem 4.2 The interpolation matrix polynomial of the Lagrange interpolation problem 
is given by 



P(x) = J2F(xi) (0 ... v iA . . . Vi, mi ... 0)(V 1 + V 2 x + ... + V n x 



(V, V 2 ... V n ) = {coliXJ^J) \ 



i=l 

where 

with Vi a np x p matrix. 
Proof 

The interpolation matrix polynomial P(x) = A n _\x n ~ l + . . . + A\x + Aq is determined by 
a linear system of np 2 equations and np 2 unknows, 



P(xi) Vij = F(xi)Vij, 
In block form we get 



J 



l,...,rrii, 



k. 



(A A 1 ... A 



n-l. 



( x \ 

XJ 
XJ 2 



\ XJ 71 - 1 J 

(F(xi)v 1}1 F(xi)v 1>2 ... F(x 1 )v 1)mi F(x 2 )v 2> i ... F{x k )v k ,m k ) 
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(see proof of Theorem |4.1|) . 

This means that the coefficients of P(x) satisfy 

(A A 1 ... = 

(F(x 1 )v 1A F(x 1 )v 1)2 . . . F(xx)v l>mi F(x 2 )v 2 ,i 



■■ F{x k )v Kmk ) -{Vi V 2 ... V n ) 



or 



Aj = (F(xi)vi A F(x 1 )v 1>2 ... F(x 1 )v hmi F(x 2 )v 2>1 ... F{x k )v Kmk )V j+1 , 
for j — 0, 1, . . . , n — 1. But then we have 

P(x) = (F( Xl )v 1A . . . F( Xl )v 1)mi F(x 2 )v 2>1 . . . F{x k )v Kmk ) (y 1 + V 2 x + ... + V n x n ~ l ) 



or 



P(x) = J2 F{ Xl ) (0 ... v itl . . . v i>mi ... 0)(V 1 + V 2 x + ... + V n x n ~ l ), 



i=i 



with (mi + . . . + m,;_i) zeros before and (m i+ i + . . . + rrik) zeros after Vi >mi . ■ 

Theorem 4.3 The interpolation matrix polynomial of the Lagrange interpolation problem 
is given by 



»=i 



where 



P(x)=J2F(x t )W l (x), 

( <i \ 



WAx) 



( 



Vi,l ■ ■ ■ v i>r 



Qn(x). 



The vector wf, is the (mi + . . . + m.j_i + j)th row from V n and Q n {x) is the monic matrix 
polynomial of degree n with Jordan pair (X, J) . 

Note that the vectors Wij are left rootvectors of the matrix polynomial Q n (x) (Section 
I)- 

Proof 

The proof of this theorem consists of three parts : 

(a) JV X - V n XJ n V x = and JV, - V n XJ n Vi = V5-i, where i = 2,...,n. 

To show these relations we consider the similar standard pair (X',C\) with 



X' = (I ... 0) 



and 



Ci 



/ 





I 


.. 





\ 










I .. 




























.. 


/ 




V 


-B 


-B x 


-B 2 .. 


• —B n - 
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where Q n (x) = Ix n + B n _ 1 x n 1 + . . . + B±x + -Bo, the monic matrix polynomial 
with Jordanpair (X,J). This means that X = X'S and J = S _1 CiS with S = 
col(XJ l )™ = y (see Section 1). The matrix col{X'C\)™ =Q satisfies 



/ X' \ 

X'd 
X'd 2 



V x'd 11 - 1 J 



( i o 

/ 



o \ 





Vo o ... i) 



because 



X'd k 

= (0 ... I ... 0) with I on the (k + l)th place and k — 0, 1, . . . , n — 1, 
= (—-Bo — Bi ... — B n -i) for k = n. 

So the matrices V( are given by 

/ \ 



V' 



with I on the ith place. 



V o J 



But then we have 
/ 1 \ ( \ 



Ci 



or 











/ / \ / \ / \ / \ 



-B —Bi . . . —B n _i) 



\ / V I J 















Consequently, 



or 



V J V -B Q J V -Bq J 

dv{ - v;x'c?v; = o. 
s^dss^v; - s- l v^x'ss- l c^ss- l v{ = o 

JV r - V n XJ n Vi = 0. 



V o / 



For i — 2,3, ... ,n we have 

/o\ /oN 




7 



v o y \ij 








^— -Bo — Bi ... — -B„_i) 



/ o N 



I 

v o y 
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/ \ / \ 



I 










/ o \ 
I 





In other words 



or 



so that 



V Si-i J V ~ B *-i J V 7 

c,v: - v n x'c?v; = vu 
s^dss^vf - s-^x'ss-^ss- 1 ^ = 
jVi - v n xrvi = Vi-L 



(b) V n Q n (x) = (xl - J) (Vi + V 2 x + . . . + V n x n ~ l ) 

The monic matrix polynomial Q n (x) with Jordan pair (X, J) is given by 

Q n (x) = x n I - XJ n (V/x + V 2 x + ... + VnX 11 - 1 ). 

(see Section 1). But this implies we immediately get the required expression from 
the relations proved in (a). 



(c) Wi(x) = ^ ( v i:1 ... Vi 



( <i \ 



V 



Qn{x)- 



From the property proved in (b) we get 



(xl - jy l V n Q n (x) = (V l + V 2 x + ... + V^- 1 ). 

If we call l mi the np x np diagonal matrix with 1 on the (mi + . . . + fflj_i + j)th 
row, where j = 1,2, ... ,m,i and on all the other places and we multiply both sides 
of the equation on the left with this matrix, we get 

/ \ 



w. 



Q n (x) = I„ H (V 1 + V 2 x + ... + V n x 



n—l\ 



V o J 

Multiplication on the left of this equation with the matrix X, gives 



1 



T 



{Vi,l v i>2 . . . v itm .) 



w 



i.2 



Qn(x) 



\ Wi tmi I 



(0 ... v itl . . . v hmi ... 0)(V 1 + V 2 x + ... + VnX 11 ' 1 ), 
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such that the previous theorem gives the required expression 

t <i \ 



k j 

P(x) = F( Xi ) — ( v iA . . . v i:TJli ) 

i=l 1 



Qn{x). 



We get another representation of the interpolation matrix polynomial if we choose the 
points Xi and the corresponding vectors in a special way. This representation shall be 
very useful for the construction of quadrature formulas. 

Theorem 4.4 Let (X, J) be a Jordan pair of the orthonormal matrix polynomial P n (x), 
then the interpolation matrix polynomial of the Lagrange interpolation problem is given 
by 

k 

P(x)=J2 F (xi)Wi(x), 
1=1 

with k the number of different zeros Xi, mi the multiplicity of Xi, Vij the vectors associated 
with Xi and 



Wi(x) = ( V it i . . . V it m t ) Ki 1 



( <1 \ 



V 



Pn+1 (Xi) T D n +t Pjx) 



The rrii x dimensional matrix Ki satisfies 

( vL \ 



Ki = 



\ V i,mi ) 



K n -i(x h Xi) ( v i:1 



with 



n-l 



K n . 1 (x,y) = J2Pj(y) T P^)- 

3=0 



Note that the Jordan chains of the orthonormal matrix polynomials P n (x) have length 
1. So, if X\,X2, ■ ■ ■ ,Xk are the zeros of the orthonormal matrix polynomials P n (x) with 
multiplicity resp. m 1 ,m 2 , . . . , m^, then every x« is associated with non-zero and linearly 
independent vectors v it i, . . . , v^ mi . Furthermore J2i=i m i — n P an d the matrix col(XJ l )™ =0 
is non-singular. In other words, these are exactly the requirements which the points and 
the vectors for the Lagrange interpolation problem, have to fulfil. 
Proof 

From the identity of Christoffel-Darboux and from the fact that Vij is a rootvector of 
P n (x) associated with the zero x iy we get 



Vf d n j2Pl(Xi) T Pl(x) = _ v j Pn+M) T D n+1 P n {x) 
1=0 



(1) 



This implies that the left side of the equation, as well the right side is a vector polynomial 
of degree n — 1 . In other words the mi x p dimensional matrix 

I <i \ 

P n+1 (xi) T D n +iPjx) 



T 

\ V i,mi ) 



'n+1 1 n\ 
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is divisible by x — x $ . 
Set 



P n +i(xi) T D n+1 P n (x) 



Wi{x) = A 

V v i,mi ) 

then this matrix polynomial of degree n — 1 satisfies 

Wi(xj) v L ,j = 0, J = 1, 2, ...,m h l = l,...,i-l,i + l,...k, 

because there is no singularity for xi ^ Xi and P n (xi)vij = 0. Now we determine the 
p x rrii dimensional matrix A such that 

Wi(xi) Vij = Vij, j = 1,2,..., m^. 

From ([!]) we get the following equation : 



n-l 



1=0 



T Pn+l( x i) Dn+lPn(x 



hi 



■ V i,s ~ V iif 



T Pn+l{Xi) T Dn+lPn(%i) 



J i,s I • 



Let x approach Xi, then this becomes 



n-l 



V i,j Pl(Xi) T Pi(Xi) V ijS = -vJjPn+ifayDn+iP^frdVi^ 
1=0 

So, the matrix A satisfies 
( <i \ 



A 



Pn+l( x i) T Dn+lP n ( X i) V i,s ~ v i,si 



S 1,2,..., ?7ij, 



V V I,m, J 



or 



A 



Consequently, A satisfies 



H.sj s 1,2,..., Tfi{. 



A 



(^) 2 ,2 



( 



or 
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Bringing everything together, we get 



Wi(x) = ( v it i . . . v i:Tni ) Ki 1 



( * \ 



V V i,rrii / 



P n +i(xi) T D n+1 P n (x) 



where 



Ki 



K n -i(xi,Xi) ( . . . ^, mi ) . 



) = -jfi 



Finally we have to show that the matrix Ki is non-singular. Therefore we show that 
—Ki is a symmetric and positive definite matrix. From the fact that K n _i(xi,Xi) is a 
symmetric and positive definite matrix, we get 

( ^ \ 

-Kj = j Kn-i(xi,Xi) ( v iA 
Let w be a fflj dimensional vector, then we have 

( W-y ... W mi ) ; 

V Mm, / 



W r , 



) 



( "1 \ 



= (Wl^jl + • • • + W mi vl mi ) K n ^{Xi,Xi) (lUiUj,! + . . . + w mi v i)mi ). 

This expression is always > and only equal to if the vector 

+ . . . + w mi ^, mi 

is equal to 0. But the vectors v^i, . . . , Vi jJUi are linearly independent, so that w T (—Ki)w = 
if and only if w — 0. So — is a symmetric and positive definite matrix and taking 
the inverse of it will not be a problem. ■ 

Special case. 

If p — 1 we consider n different points x±, . . . , x n and associate a non zero number Vi with 
every point Xj. The interpolation polynomial for a function /(#) is given by 



with 



Wi(x) = -v. 



i=l 

Pn+l(^i)^n+lPn(^) 



-1 

ViK n _\(Xi^ Xi^Vi X Xi 

~Pn+l( x i)k n ,n Pn( x ) 
kn+l,n+lK n —\{Xi, Xi) X Xi 

Pn(x) 
(X - Xi)p' n (Xi)' 
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Example. 

Consider the orthonormal matrix polynomials on the interval [—1,1] with respect to the 
weight matrix function 



W(x) 



i(l_a;2)-V2 o 

o i(i-* 2 ) 1/2 



The recursion coefficients are given by : 

E n = for n > 0, 

£>i = ( ^ i ) and °n = ( ^ i ) for n > 1. 
In this case the orthonormal matrix polynomials are 



V2x \ 

o 2x r 



4x 2 -lJ' 8x 3 -4xJ- 

In general, we have 

- ( * T u! (x) ) ■ 

Consider p2(x), then the zeros and the corresponding rootvectors are given by 

1 -11-1 



n / 4 \ o\ f o 

This leads to 



^ = ' U2 = ' ^ = 3 ' " 4 ~ V -2 



(;) — » 



and 



^00 

In a similar manner we find 



K 2 = -32, W 2 (x) = ( ^ ] ° 



or 



X 3 = -18, W 3 (x) = 







(x + l) 
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K t = -8, W t ( X ) = ( I _ {x °_ , } ) . 

Let F(a;) be given by 

F , v / 2a; + 5 Qx \ 
F W=[ 7 Ax-3 J' 

then the interpolating matrix polynomial is given by 

/ 2a; + 5 Qx \ . , 
7 4a;-3 = ^ 



Let be given by 



x 2 + 1 6a; 



x 7a; + 1 5a; 2 - 1 J ' 
then the interpolating matrix polynomial is given by 

P(x) = ( x 3 f ) W x {x) + ( _J § + ! " 3 |^ ) W 2 (x) 
/ 5 3 \ TTr , x / 4 -3 x 



+ (J i J W 3 (x) + ^ j W 4 (a;) 

/ | 6a; ^ 
7a; + 1 i 



4 



5 Gaussian quadrature 

Let W(x) be a matrix weight function defined on the interval [a, b]. Then we are going to 
approximate the integral of matrix functions by means of a sum of the form 



b k 

( F(x) W(x) G(x) T dx~Y. F ( x *) A * G ( x if 
Ja i=i 



where A,- G R pxp . 



5.1 Quadrature formulas 

The numerical integration problem exists in finding points Xi and matrices Aj, where 
% — 1, . . . , k, so that the formula will have the highest possible degree of precision. This 
degree of precision of a quadrature formula is equal to m(e N) if 

f x l W(x)dx-Y J A A i = for Z = 0,1,..., m, 

Ja i=i 

^0 for Z = m + 1. 
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Note that this agrees with the accuracy of the quadrature formula 

rb k 



f F(x) W{x) G(x) T dx~f^ F ( x i) A * G (^) T > 
Ja i=i 



for matrix polynomials which satisfy degF(x) + deg G(x) < m. 
It will be convient to choose 



K = ( 



Ai 



( <i \ 



T 

\ V i,rrn ) 



with Vi t i, . . . , Vi >mi linearly independent, non-zero vectors and Yn=i m % — n V- The numeri- 
cal integration problem exists than in determining Xi, vectors Vij and matrices A4, where 
j = 1, . . . , m,i and % — 1, . . . , k, so that the formula will have the highest possible degree of 
precision. We will show now that 2n — 1 is the highest degree of precision for a quadrature 
formula of the postulated form. 

Theorem 5.1 A quadrature formula of the form 

f b F(x) W{x) G(xf dx~J2 F{xi) A* G{ Xl ) T 
Ja i=i 

cannot be exact for all matrix polynomials F(x) and G(x) satisfying 

degF{x) + degG{x) < In. 

Proof 

Let F[x) = G(x) = Q n (x), where Q n (x) is the monic matrix polynomial of degree n which 
satisfies Q n (%i)vij = for j = 1, 2, . . . , and i — 1, 2, . . . , k. In this case we have 

f F{x)W{x)G{x) T dx= f b Q n (x) W(x)Q n (xfdx > 0. 

Ja Ja 



On the other hand we have 



J2F(xi)AiG(xi) T = J2Qn(xi) ( v i:1 ... v ijmi ) Ai 



( <i \ 



i=i 



i=l 



Qn(Xi) = 0, 



\ 4 / 



so that the quadrature formula cannot be exact in this case. 



In the following properties we will show that it is possible to construct a quadrature 
formula with degree of precision 2n — 1. 

Theorem 5.2 Let (X, J) be a Jordan pair of the orthonormal matrix polynomial P n {x) 
on the interval [a,b] and with respect to the weight matrix function W(x). Then we have 



f b F(x) W(x) G(x) T dx~Y. F ( x i) A * G ( x i? 
Ja i=i 



(2) 
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where k is the number of different zeros %i, mj is the multiplicity of X{, Vij are the vectors 
associated with x i; 

( <i \ 

Ai = ( • • • v i>mi ) L~ l 



V V i,mi j 



and 



v 4 J 



K n _i(xi,Xi) ( v iA 



This quadrature formula is exact for matrix polynomials F(x) and G(x) which satisfy 

deg F(x) + deg G(x) < 2n - 1. 

Proof 

If we replace F(x) and G(x) by their interpolating matrix polynomial from the Lagrange 
interpolation problem (see Theorem |4.4| ), then we get 

njj k k 

/ F{x) W(x) G(x) T dx~J2Yl F ( x i) A m G ( x jf> 

Ja i=l 7=1 



with 



( <i \ 



A; 



'•j 



( • • • f«,m, 

J a 



K- 1 



b P (x) 
— — W[x / 



PJx) T 



dx D n+1 



( 



X V 3,l 



) Kf 1 



( v l \ 



From the identity of Christoffel-Darboux we get 



( <i \ 



v i, / 



( <i \ 



K n —i(x, Xi 



P n+l { Xl ) T D n+l P n {x) 



\ %m J 



so that 

Aij = ( v iA 



( <i \ 



J i,rrn 



Kr 1 



V v i,m,i J 



rp 

K n -i[x, Xi) W(x) K n -i(x, Xj) dx 



( »h \ 



X Vi,l 



Vi,l . ■ ■ fi,r 



AV 1 



V V lm, J 
I <1 \ 



K 



-i 



t * \ 



v 4 j 



K n -i{xj,Xi) ( v jA 



u 3,m i j Kj 



-1 
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But for % ^ j we have 



Pn{%i) T D n+ iP n+ i(Xj) — P n+ i(Xi) T D n+ iP n (Xj) 



Xj X 2 



V 3,t 



o, 



so that Aj j = for i 7^ j. For i = j we get 



Aj = A M = - ( u^x . . . v i:Tni Ki 1 KiKi 1 



( <i \ 



where 



—Ki 



( v iA 



V V i,rrii ) 



K n _i(x h Xi) ( 



which was to be proved. 

Now we show that this quadrature formula is exact for matrix polynomials F(x) and G(x) 
satisfying 

degF(x) + degG(x) < 2n - 1. 

If F(x) and Cr(x) are matrix polynomials with degree < n—1, then the quadrature formula 
will be exact because the interpolation matrix polynomials are exactly F(x) and G(x). 
Let F(x) be a matrix poynomial of degree n + I and a matrix polynomial of degree 
n — I — 1, where Z = 0, 1, . . . , n — 1. Then Gr(x) satisfies the equation 

G{x) = Y J G{x i )W l {x), 
i=i 

because deg < n — 1. On the other hand F(x) is approximated by a matrix polyno- 
mial of degree n — 1 : 



~^F(x,)W,(a;) = P( 



i=l 



This implies that F(x)—P(x) is a matrix polynomial of degree n+/, with zeros X\, x 2 , ■ ■ ■ , x^ 
and rootvectors v^i, . . . , i>i, mi , . . . , i^i, • • • , v k,m k - But then we have 

F{x)-P{x) = R{x)P n (x), 

with R(x) a matrix polynomial of degree I, (see Section 1). So we get 

[ b F(x)W(x)G(x) T dx = [ b P(x)W(x)G(x) T dx + f R(x)P n (x)W(x)G(x) T dx 

J a J a J a 

= j2P{x t )KG{ Xl ) T + tij^Rix'] P n (x)W(x)G(x) T dx 

i=l Ja \i=0 ) 
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k I b 

Y J P{x i )K i G{x l ) T + Y,Ri [ x i P n {x)W{x)G(xfdx 

i=l i=0 Ja 

J2P{xi)KG{xi) T + Y,Ri t Pn{x)W{x)x i G{x) T dx 



The last integral will vanish because degx l G(x) < I + n — I — 1 = n — 1 and P n (x) is 
orthogonal to every matrix polynomial of lower degree. Further we have 



P(Xi) A; = ( 



F(xi) ( 



• • • ui,mi 



LT 1 



LT 1 



( <l \ 

\ V lmi J 
( <1 \ 



V V i,rrii ) 



so that 



b k 

/ F(x) W(z) G(x) T dx = J2 F{xi) K G(xi) rj 
Ja i=i 



The case where deg G(x) > n is treated in a similar way, but we can also take the transpose 
of the above mentioned equation. ■ 

We showed that using the zeros and the rootvectors of the orthonormal matrix polynomials 
on the interval [a, b] with respect to the weight matrix function W(x) leads to a quadrature 
formula with degree of precision 2n — 1. 



Note that Lj is a symmetric and positive definite matrix (see proof of Theorem 4~4] ) and 

so Aj is a symmetric and positive semidefinite matrix. 

Example. 

Consider the orthonormal matrix polynomials on the interval [—1,1] with respect to the 
weight matrix function 



W(x) = 
We already mentioned that 

P2(X) 

with zeros 



±a-* 2 r 1/2 





1 — x 



2\V2 



2^2x 2 -y/2 







Ax 2 



1 



Xl 



V2' 



X2 



V2' 



;r.i 



and rootvectors 



1 




^2 



4 




1>3 





3 



1> 4 





-2 
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The quadrature coefficients can be computed as follows 



-l 



i (-) r.U i)) 



1/2 

o o 



In a similar way we find 



1/2 \ / \ / 

J A3 ~ 1 1/2 ) A4 " { 1/2 



Suppose we have the following matrix polynomials 

F(a;) = ( 7x + 1 5* 2 -l J and G ^ = { 7 4x-3 

then 



3 



-3^ \ . / 5- V2 7 



7 + l I j A2 I -3V2 -2V2-3 



2 



/ 33/2 33/2 \ 

^ 12 25/4 ) 



5.2 Convergence 

The Gaussian quadrature formula converges to the exact value of the matrix integral of 
the functions F(x) and G(x), even without imposing severe conditions on those functions. 

Theorem 5.3 Let F(x) and G(x) be continuous matrix functions on the finite interval 
[a, b], then we have 

f F(x) W(x) G{x) T dx = lim ^F(^ n) ) A 4 (n) G^f , 

J a n ~* co i=1 

with x\ n ^ the zeros of the orthonormal matrix polynomial P n (x) on the interval [a, b] with 
respect to the weight matrix function W(x) and with the corresponding quadrature 
coefficients. 

In the proof of this theorem we will use the following proposition of matrix norms: 
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Proposition 5.4 If Ai,A 2} . . ., are symmetric positive semidefinite p x p-matrices, 
then they satisfy 

E\\a\\ 2 <p X> 

with 

\\A\\ 2 = ^p(A*A) = p(A) 
and p(A) is the spectral radius of A. 
Proof 

Suppose a® > ... > > are the eigenvalues of A4 and 71 > • • • > 7 P > are those 
of Yjl=i Ai- Then we have 



E 

i=l 



E«?<EE«f = EM^)=fr£^ 



..w 



* 1 1 2 



E^i < P7i = P 



i=l j'=l i=l 
k 



i=i 



i=i 



If the interval [a, b] is finite we have another useful inequality, namely 

|| f F(x)dS(x)\\ < ! \\F{x)\\d8{x), 

J a J a 

PL p. 78]). 



sec 



Proof of Theorem 5.3 



We will distinguish 3 cases : 
(1) G{x) =1. 

Since the matrix function F(x) is continuous, every Fij(x) will be continuous. This 
means that every element can be approximated arbitrarily close by a polynomial 
and thus 

\\F{x)-Q m {x)\\ 2 <e. 



So we have 



f F(x) W(x) rfx-E F (4 n) ) AS n) I 
Ja i=i 



< 



F(x)W{x)dx- Q m (x)W(x)dx\\ 2 



+ II f Q m {x)W{x)dx-j2Q m {x i T ) )^r ) h 
Ja i=l 



[n)\ A ( n )| 



i=l 



i=l 



But the first term of this sum satisfies 



|| [ b F(x)W(x)dx- f Q m (x)W(x)dx\\ 2 < f \\F{x)-Q m {x)\\ 2 \\W(x)\\ 2 dx 

J a J a 

< € [ b \\W(x)\\ 2 dx. 



26 



The second term 

rb 



( Q m (x) W{x) dx-Y, Q m (4 n) ) AS n) || 2 = 0, 
Ja i=l 



for m < In — 1 and the third term can be bounded as follows : 

11 Y.q m {xf) a<"> - tn4 n) ) M n) h < e ng m (^ (n) ) - Fixt^h \\^\ 

i=l i=l i=l 

k 



< ^EiiAr'ib. 



i=i 

So we have 

r 6 fc , , , , / r b k 



||/ F{x)W{x)dx-Y J F{4 l) )^ l) h<e( I \\W{x)\\ 2 dx + Y,\\^ 
Ja i=\ \ Ja i=i 

But from the previous proposition we get 

El|A! n) || 2 <p|lEA! n) || 2 =p|| f b W(x)dx\\ 2 <p f \\W(x)\\ 2 dx 

■1=1 i=i Ja Ja 



(«)n 2 



and thus we have 

r b 



f F(x)W(x)dx-^2F(x\ n) )Al n) \\ 2 <e(l+p) I \\W{x)\\ 2 dx. 

J a i=l 



But this means that 

rb 



f F(x)W(x)dx= lim E^S n) )Aj n) . 



(2) is a matrix polynomial. 

Suppose G(x) = G;x' + + . . . + Gix + G , then we have 



/ F(x)W(x)G(x) T dx = E(/ G' 

a i=0 y a J 

= E (j^x i F(x)W(x)dx\ G) 



3=0 

= EilmE(^ n) )^ri n) )AS n) Gj 
i=o i=l 

= ilmE^^Al^E^^J 
i=l i=o 

= iimE^ri n) )AS n) G(^f. 
i=i 

Which implies that the theorem holds in the case that G(x) is a matrix polynomial. 
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(3) G(x) is a continuous matrix function on the interval [a, b]. 

In this case, both functions F(x) and G(x) can be approximated arbitrarily closely 
by a matrix polynomial : 



This gives 



+ 



+ 



+ 



+ 



\\F(x)-Q m (x)h<e and \\G(x) - P,(a;)|| 2 < e. 

f F{x) W(x) G(x) T dx-f: F{x [ t ] ) Al n) G(x[ n) f\\ 2 
Ja i=i 

f F{x) W{x) Gix) T dx - f Q m ix) W(x) Gix) T dx 

la Ja 

- f ' F{x)W{x)P l {x) T dx+ f ' Q m ix)Wix)P l ix) T dx\\ 2 

Ja Ja 

fQ m ix) W(x) Gixfdx -j^QMf) A™ G(x<r>f\\ 2 

Ja 1=1 

f b Fix) W(x) Ptixf dx-^2 Fix^) A, (n) Ptix^fh 
Ja i=l 

fQ m ix) Wix) Vixfdx -j2Q m (x\ n) ) A, (n) Pi(x[ n) fh 

Ja i=l 



Y.Qmixf) ^Gixf) +Y,n4 n) ) aS" } Piix\ n) ) 



1=1 

k 



1=1 

J n )\ a(«) r>(A n )\ T EY^( n h A (") r<(„( n )^ T 



1=1 



1=1 



From the previous case we get that the second and the third term are arbitrarily 
small (< e) for n sufficiently large. The fourth term is equal to if m + I < 2n — 1 
and the other two terms can be bounded as follows, 

|| ["[Fix) - Q m (x)} Wix) [Gix) - P^fdxh 

J a 

< t ||F(a;)-g m (a;)|| 2 ||W(a;)|| 2 ||G(a;)-P I (x)|| 2 tfa; 



< 



e 2 f \\W{x)\\ 2 dx 

J a 



and 



|| - Qm(x\ n) )\ Ar [G{xf } ) - Pl(xF)]* \\ 2 



Hm a H \n(J n h 



i=i 



< E H^ (n) ) - QmixDh \\kT% \\g{xT } ) - Pi(x\ n, )h 



Wml ii a Hi 



i=l 



i=i 
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Hence for n large enough 

|| f F{x) W(x) G(x) T dx - J>(z< n) ) A, (n) G(xS n) ) T || 2 
■ /a i=i 

< 2e + £ 2 (^V(*)M* + E||Af ) || 2 ) 

< 2e + e 2 (l+p) [ b \\W(x)\\ 2 dx, 

J a 

which proves the convergence of the Gaussian quadrature formula. ■ 
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